Interpretation of Hund's multiplicity rule for the carbon atom 
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Abstract 

Hund's multiplicity rule is investigated for the carbon atom using quantum Monte Carlo methods. 
Our calculations give an accurate account of electronic correlation and obey the virial theorem to 
high accuracy. This allows us to obtain accurate values for each of the energy terms and therefore 
to give a convincing explanation of the mechanism by which Hund's rule operates in carbon. We 
find that the energy gain in the triplet with respect to the singlet state is due to the greater 
electron-nucleus attraction in the higher spin state, in accordance with Hartree-Fock calculations 
and studies including correlation. The method used here can easily be extended to heavier atoms. 



I. INTRODUCTION 



From an analysis of atomic spectra, Hund found that the electronic configuration of the 
lowest energy is the one with the highest spin multiplicity S. Hund's multiplicity rule 3] 
explains almost all the configurations of the ground states of atoms, ions, and molecules 
and even their low-lying excited states in most cases. Even simple models such as frozen 
orbital Hartree-Fock (HF) theory generally adhere to the rule, although the mechanism by 
which Hund's rule is obeyed depends on the accuracy of the theoretical description. A widely 
accepted mechanism given originally by Slater attributes it to a reduction in the electron- 
electron repulsion energy V^g in the highest S state. This mechanism is in fact incorrect as it 
is based on the incorrect assumption of choosing the same orbitals for states with different S 
(frozen orbitals). This assumption leaves the kinetic energy T unchanged for different spin 
states, which cannot be correct as each eigenstate must satisfy the virial theorem. The virial 
theorem states that 2T+V = 0, where V is the total potential energy, which shows that if the 
potential energy is decreased the kinetic energy must increase. If the orbitals are allowed to 
relax self-consistently they are found to depend on 5* states. 

The correct explanation of Hund's multiplicity rule is in fact that the higher spin state 
reduces the electron-nucleus interaction energy Ven- Interestingly it turns out that states 
with larger S have higher values of Vee, in direct contradiction to the Slater mechanism. 

Davidson 4] confirmed that Hund's multiplicity rule derives from a reduction in the 
electron-nucleus energy within the self-consistent HF approximation for the low-lying excited 
states of the helium isoelectronic series. The self-consistent HF method has been applied 
to explain Hund's multiplicity rule for all atoms (with the exception of Zr (Z=40) for 



IJ], giving the same 



which HF theory predicts an incorrect ground state) 

stabilizing mechanism. Boyd0, |l^ considered the same systems as Davidson but within a 
variational approach. He explained why larger values of 5* give lower energies as follows. 
Larger 5* implies more parallel spin pairs and hence a larger exclusion holes around each 
electron. Consequently, electrons can avoid each other more effectively and move closer 
to the nucleus, which reduces V^n- We refer to this as the less screening mechanism. We 
should expect that when a mechanism lowers Ven it should also increase Vee, because the 
electrons will on average be closer to one another. The less screening mechanism is therefore 
in contradiction to Slater's argument which implies that Vee will be lower in the triplet state. 
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Configuration interaction (CI) calculations have shown that the same mechanism operates 



in the low-lying excited states o. 
hydrides, LiH, BH, NH, and FH. 



small molecules such as CH2, H2CO, and the first-row 
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The origin of Hund's multiplicity rule for heavier atoms has mostly been examined at 
the HF level. It is interesting to study the origin of Hund's multiplicity rule in heavier 
atoms beyond the HF level. However, the conventional approaches such as CI are difficult 
to apply for the following reasons: (1) The wave function must be constructed carefully from 
the configuration state functions (CSFs) giving the largest contributions to the correlation 
energy, which becomes a cumbersome task for heavier atoms. (2) A huge number of CSFs are 
required to give a sufficiently good description of the correlation, particularly the electron- 
electron cusps. In this paper we report a promising approach to this problem applied to the 
carbon atom as a prototype using quantum Monte Carlo (QMC) approaches such as the 
variational Monte Carlo (VMC) and more accurate diffusion Monte Carlo (DMC) methods. 

II. DESCRIPTION OF THE CALCULATIONS 

In VMC, expectation values are calculated using an approximate many-body wave func- 
tion, the integrals being performed by a Monte Carlo method. The approximate wave 
function normally contains a number of variable parameters, whose values are obtained by 
an optimization procedure. 

DMC is a stochastic projector method for solving the imaginary-time many-body 
Schrodinger equation. Exact imaginary-time evolution would in principle give the ex- 
act ground-state energy, provided the initial state has a non-zero overlap with the exact 
ground state. However, the stochastic evolution is not exact and the symmetry of the final 
state would not be fermionic, but bosonic. To maintain the correct fermionic symmetry the 
fixed-node approximation jly, |l2| is used in which the nodal surface of the wave function 
is constrained to equal that of an approximate guiding wave function. If the nodal surface 
of the guiding wave function is accurate, fixed-node DMC will give a very accurate energy. 
The fixed-node DMC energy is less than or equal to the variational energy calculated with 
the guiding wave function and greater than or equal to the exact energy. 

We used a single-determinant Slater- Jastrow guiding wave function, jisl The single- 
particle orbitals were obtained from HF calculations using the GAUSSIAN98 code and 
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a 6-311G++(3df,2pd) basis set. At a nucleus the exact HF orbitals have cusps such that the 
divergence in the potential energy is cancelled by an equal and opposite divergence in the 
kinetic energy. Orbitals expanded in a Gaussians basis set cannot have cusps and the sum of 
the contributions from the kinetic and electron-nucleus energies diverges at the nucleus. In 
practice one finds that this results in wild oscillations in the local energy when an electron 
approaches the nucleus, which can lead to numerical instabilities in DMC calculations. To 
solve this problem we make small corrections to the single particle orbitals close to the nuclei 
which impose the correct cusp behavior. 

Electron correlation is introduced by multiplying the determinant by a Jastrow function 
which depends explicitly on the inter-electronic distances and obeys the electron-electron 
cusp conditions. |22] Our Jastrow factors take the form 



i>j 



A 



expi-- 



exp 



ij 



where the indices i and j denote electrons and / denotes ions. F is chosen so that the cusp 



conditions 



22j are obeyed, i.e., F-^ 



TT 



2A and Fi 



n 



A, and Si and S2 are cusp-less. 



5*1 and S2 are expressed as polynomial expansions in the inter-particle distances. We used 
a total of 12 variable parameters in the Jastrow factor whose optimal values were obtained 



by minimizing the variance of the VMC energy. The wave functions obtained by this 

procedure were used as the guiding wave functions for our DMC calculations. We used a 
timestsep of 0.001 a.u. within the DMC calculations, which gives ve ry s mall timestep errors. 
All the QMC calculations were performed using the CASINO code. |21 1 

Arguments about the mechanism by which Hund's rule operates centre on the values of 
the individual energy terms T, V^e and Ven- As illustrated by Slater's early treatment, results 
which do not obey the virial theorem to high accuracy can give a spurious interpretation of 
the rule. It is therefore necessary to calculate the energy terms using an accurate method, 
and in the case of a stochastic method such as DMC, to evaluate them to a high statistical 
precision. Because there is no "zero-variance property" for the individual energy terms it is 
necessary to perform long Monte Carlo runs, and in our calculations we have accumulated 
the results over 200,000 steps. 

The DMC method generates the mixed distribution, \&t,$05 where and $0 are the 
guiding (VMC) and DMC wave functions, respectively. Evaluating the expectation value of 
the Hamiltonian H (or any operator which commutes with it) with the mixed distribution 



gives an unbiased estimate of the energy, i.e., 

'6) -™^-^^^H^> . P) 



($0 




($0 





p 



where the subscript m denotes the mixed estimator and p denotes the pure estimator. 
However, for expectation values of operators which do not commute with the Hamiltonian 
the mixed estimate is different from, and generally less accurate than, the pure estimate. An 
approximation to the pure estimate can be obtained by combining VMC and DMC results 
in an extrapolated estimation, j^] The linear extrapolated estimator is 

{d)i = 2{d)m - (O). = (0)p + 0(A2), (3) 

where A = $0 ~ ^t)- Other extrapolated estimators such as the squared extrapolation, 

(6). = ^ = (6), + 0(A2), (4) 

can also be used, which have different error terms. In this study we found that the linear 
and squared extrapolations gave almost identical results, and we will quote results only for 
the linear one. 



III. RESULTS AND DISCUSSION 

In Table I we give VMC and DMC expectation values for the different each energy terms 
and for the virial ratio for the lowest energy triplet and sing let states of the carbon atom. 



For comparison, the HF and CI values of Pauncz et al. 6, 



14| are also given. The exact non- 



relativistic energy of the ground state of the carbon atom has been estimated to be -37.8450 



a.u. 3 



We can use this value and the Hartree-Fock energy of -37.689 a.u. to compute the 



percentages of the correlation energy recovered. The CI calculations of Pauncz et al. 

n 

recovered 45.6 % of the correlation energy, while a more modern CI calculation [28] (which 
did not report the individual energy terms) recovered 61.5 %. Our DMC calculations recover 
88.5 % of the correlation energy. 

For each of the VMC, DMC and extrapolated DMC calculations we obtain a virial ratio 
within statistical error bars of its correct value of —V/T = 2. It is straightforward to show 
from the virial theorem that the exact kinetic energy should be larger than the HF one in 
an isolated system. This inequality is obeyed by each of our calculations but is disobeyed by 



the CI calculations of Lemberger and Pauncz,[14i] whose calculations also show a significant 
error in the virial ratio. This indicates that the virial ratio provides a sensitive test of 
whether the relationships between the different energy components are correct. 

Comparing the DMC and extrapolated DMC as results with the HF ones for a particular 
state we find 

V < V^^ 

yen ^ I' en 

T > T"^. (5) 

The introduction of correlation reduces both the electron-electron and electron-nuclear po- 
tential energies at the expense of increasing the kinetic energy. This increase can be un- 
derstood as a consequence of the uncertainty principle; the increased correlation allows the 
charge density to contract towards the nucleus, which lowers V^n- However, the more local- 
ized charge density leads to an increase in the kinetic energy. Such a mechanism is not well 
described by the DCI (doubly excited CI) wave functions employed in Pauncz's calculations 
because the single excitation terms which are excluded contribute to the modification of the 
charge density. The spurious decrease in the kinetic energy seen in Pauncz's results might 
be attributed to this. In the QMC approach, on the other hand, such a modification of 



the charge density can be described by the one-body term in the Jastrow function [zyj even 
within single determinant calculations. 

Comparing the DMC and extrapolated DMC results for the singlet and triplet states we 
find 



s=i ^ yS=o 

S=l ^ T/'S'=0 

m en 

S=l ^ rp5=0 



Both the electron-electron repulsion and kinetic energy are lower in the singlet state, and 
the stabilization of the triplet occurs via the larger reduction in the electron-nucleus energy, 
corresponding to a contraction of the electron density towards the nucleus. This contraction 
is promoted by the presence of more spin-parallel pairs in the triplet which enables the 
contraction to occur without a large increase in the electron-electron interaction energy. 
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The higher values of Vee and T in the triplet state are a consequence of the contraction 
depending on the spin state, which cannot be described within Slater's frozen orbital model. 

Our results, which include a high-level description of correlation and satisfy the virial 
theorem to high accuracy, support Boyd's less screening interpretation of Hund's multiplicity 
rule for the carbon atom. It would not be difficult to extend this study to heavier atoms as, 
for example, all-electron DMC calculations of Xe (Z=54) have recently been performed, j^] 
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TABLE I: Energies for the triplet and singlet states of the carbon atom. All energies are given in 
hartrees. Standard errors in the QMC results are indicated in parentheses. 



1\ /r^-i- V. ^ /I 
Metnoa 


btate 


Energy 


Vcc 


T/ 

Vcn 


1 


Virial ratio 




Triplet 


-37.689 


12.760 


-88.137 


37.689 


2.000 




Singlet 


-37.631 


12.728 


-87.992 


37.632 


2.000 




Triplet 


-37.760 


12.563 


-87.956 


37.635 


2.003 




Singlet 


-37.703 


12.535 


-87.821 


37.583 


2.003 


VMC 


Triplet 


-37.7717(1) 


12.629(1) 


-88.174(13) 


37.773(12) 


2.000(1) 




Singlet 


-37.6896(1) 


12.553(1) 


-87.932(13) 


37.689(13) 


2.000(1) 


DMC 


Triplet 


-37.8267(4) 


12.583(3) 


-88.203(13) 


37.793(12) 


2.001(1) 


(mixed estimator) 


Singlet 


-37.7623(6) 


12.509(5) 


-87.997(22) 


37.726(19) 


2.001(1) 


DMC 


Triplet 




12.538(6) 


-88.231(29) 


37.811(26) 


2.002(2) 


(linear estimator) 


Singlet 




12.464(9) 


-88.062(45) 


37.763(41) 


2.002(2) 



Reference [6j ^Reference fl4l 
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